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Abstract. The simplest unified extension of tlie Minimal Supersymmet- 
ric Standard Model with bi-linear R-Parity violation provides a predictive 
scheme for neutrino masses which can account for the observed atmospheric 
and solar neutrino anomalies. Despite the smallness of neutrino masses R- 
parity violation is observable at present and future high-energy colliders, 
providing an unambiguous cross-check of the model. 



1. Introduction 

The recent announcement of high statistics atmospheric neutrino data by 
the SuperKamiokande collaboration [1] has confirmed the deficit of muon 
neutrinos, especially at small zenith angles, opening a new era in neutrino 
physics. Although there may be alternative solutions of the atmospheric 
neutrino anomaly [2] it is fair to say that the simplest interpretation of 
the data is in terms of to z/,- flavour oscillations with maximal mixing. 
This excludes a large mixing among I'r and ly^ [1], in agreement also with 
the Chooz reactor data. On the other hand the persistent disagreement 
between solar neutrino data and theoretical expectations [3] has been a 
long-standing problem in physics. Recent solar neutrino data [4] are con- 
sistent with both vacuum oscillations and MSW conversions. In the latter 
case one can have either the large or the small mixing angle solutions, with 
a slight trend towards the latter [5]. 

Many attempts have appeared in the literature to explain the data. 
Here we review recent results [6] obtained in a model [7] which is a simple 
extension of the MSSM with with bilinear R-parity violation (BRPV). This 
model, despite being a minimal extension of the MSSM, can explain the 
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solar and atmospheric neutrino data. Its most attractive feature is that 
it gives definite predictions for accelerator physics for the same range of 
parameters that explain the neutrino data. 

2. Broken R— parity 

In the past most discussions of supersymmetric (SUSY) phenomenology 
assumed R-parity (Rp) conservation where, 

Rp = (_1)2J+3S+L 

This implies that SUSY particles are pair produced, every SUSY particle 
decays into another SUSY particle and that there is a LSP that it is sta- 
ble. But this is just an ad hoc assumption without a deep justification. In 
this talk we will review how Rp can be broken, either spontaneously or 
explicitly, and discuss the most important features of these models [8] . 

2.1. SPONTANEOUSLY BROKEN R-PARITY 

2.1.1. The Original Proposal 

In the original proposal [9] the content was just the MSSM and the breaking 
was induced by 

{i>r) = VL (2) 

The problem with this model was that the Majoron J coupled to with 
gauge strength and therefore the decay — > p^J contributed to the in- 
visible Z width the equivalent of half a (light) neutrino family. After LEP 
I this was excluded. 

2.1.2. A Viable Model for SB RP 

The way to avoid the previous difficulty is to enlarge the model and make 
J mostly out of isosinglets. This was proposed by Masiero and Valle [10]. 
The content is the MSSM phis a few Isosinglet Superfields that carry lepton 
number. The model is defined by the superpotential [10], 

W = huu''QHu + hdd^QHd + hee''LHd 

where the lepton number assignments are shown in Table (1). The sponta- 
neous breaking of R parity and lepton number is driven by [10] 



VR = {i^Rr) VS = (Sr) VL = (j>r) 



(4) 
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TABLE 1. Lepton number assignments. 

Field L e" S others 

Lepton #1-1-11 



The electroweak breaking and fermion masses arise from 

{Hu) = Vu {Hd) = vd (5) 

with v"^ = + v"^ fixed by the W mass. The Majoron is given by the 
imaginary part of 

1!^ V V V 

;{VuHu-VdHd) + —l7r - + —Sr (6) 



where V = y v'j^ ~^ "^s- Since the Majoron is mainly an SU (2) C/(l) singlet 
it does not contribute to the invisible decay width. 

2.1.3. Some Results on SBRP 

The SBRP model has been extensively studied. The implications for accel- 
erator and non-accelerator physics have been presented before and we will 
not discuss them here [8]. In this talk we will only review the results for 
neutrinos. Neutrinos are massless at Lagrangian level but get mass from the 
mixing with ncutralinos[ll]. Neutrinos mix and the mixing is related to the 
the coupling matrix h^^. . This matrix has to be non diagonal in generation 
space to allow 

Vr^V^ + J {7) 

and therefore evading [11] the Critical Density Argvment against u's in the 
MeV range. In the SM BBN arguments [12] rule out ly-r masses in the range 

0.5 MeV < < 35MeV (8) 

We have shown [13] that SBRP models can evade that constraint due to 
new annihilation channels 

Ut-Vt — > J J (9) 



2.2. EXPLICITLY BROKEN R-PARITY 

The most general superpotential W with the particle content of the MSSM 
is given by [7] 
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W = Wmssm + Wfi 



(10) 



where 

Wmssm = e^b [h^^QtuM + h%Q\DM + h%L%H2 - fiH^Ht] (11) 
and 

Wfi = Sab [Xijk^L'^jRk + KjA^Ql + KjADjUk + Sab eiLtHt] (12) 

where i,j = 1,2,3 are generation indices, a,b = 1,2 are SU{2) indices. To 
these we also have to add the soft supersymmetry breaking terms [6]. 

3. Bilinear R-Parity Violation (BRPV) 

3.1. THE MODEL 

The superpotential W is given by 

W=eab[h'i}QtUjHt+ h%Q\b^H2+ h%L\RjH2- iimHt+ eiZfHt] (13) 
while the set of soft supersymmetry breaking terms are 



kMiXiXi + Sab {A'^QtUjH: 



+A%Q\D^H2 + A%L\R^H2 - B^iH^H'^ + B,e,Lf + h.c 



(14) 



The bihnear R-parity violating term cannot be eliminated by supcrfield 
redefinition. The reason is [14] that the bottom Yukawa coupling, usually 
neglected, plays a crucial role in splitting the soft-breaking parameters B 
and Bi as well as the scalar masses rnjj^ and M£, assumed to be equal at 
the unification scale. 

The BRPV model is a 1(3) parameter (s) generalization of the MSSM. It can 
be thought as an effective model showing the more important features of 
the SBRP-model at the weak scale. The mass matrices, charged and neutral 
currents, are similar to the SBRP-model if we identify 



e = VRhi, 

The model has the MSSM as a limit when e,- — 0. 



(15) 
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3.2. RADIATIVE BREAKING 

At (5 = Mqijt we assume the standard minimal supergravity unifications 
assumptions, 



At = A,, = Ar = A,B = B2 = A-l, 

ml, = ml^ =Ml = Ml = Ml = Mlj = Ml = ml, 

Ms = M2 = Ml = Ml/2 



(16) 



In order to determine the values of the Yukawa couplings and of the soft 
breaking scalar masses at low energies we first run the RGE's from the 
unification scale Mqut ~ 10^^ GeV down to the weak scale. For details 
see [6, 7]. 



4. Tree Level Neutrino Masses and Mixings 
4.1. NEUTRAL FERMION MASS MATRIX 

In the basis ilP'^ = (— iA', —iX^, H^, H^, Vg, f^, Vt-) the neutral fermions mass 
terms in the Lagrangian are given by 
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where the neutralino/neutrino mass matrix is 

Mjv = 







with 



Ml 


-\9'vd 



M2 

\gvd 
-\gvu 



m 



\g'»d -\gvu 
-/i 



(17) 



(18) 







ra 



0.2 
L as 



(19) 



where = {—\g'vi^\gvi^Q^ei). This neutralino/neutrino mass matrix is 
diagonalized by 



H* MnM ^ = diag (jTij^o , m^o , m^o , m^o ,mi,^,mi,^,mi,^) 



(20) 



4.2. APPROXIMATE DIAGONALIZATION OF MASS MATRICES 
If the l^p parameters are small, then 

^ = m • M-^l => y^ij < 1 



(21) 
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one can find an approximate solution for the mixing matrix M. Explicit 
solutions can be found in Ref. [6]. In leading order in ^ the mixing matrix 
is given by, 



N* 




1 - m 

1 - 



(22) 



The second matrix above block-diagonalize Mjv approximately to the form 
diag(me//, A^^o) 



me// = -m ■ Mjm^ = . \ AgA^ Af, A^A^ 



= Tr{meff) 



idet{Myo] 



(23) 



The submatrices N and V^, in Eq. (22) diagonalize M^o and nieff 

N*M^oN^ = diag{m^o) ; VjnieffV^ = diag{0,0,m„), (24) 

where 



(25) 



5. One Loop Neutrino Masses and Mixings 

5.1. DEFINITION 

The Self-Energy for the neutralino/neutrino is 




^=i{]^ [Pi^fj + Pr^S\ - [Pl^^J + Pr^S\ } (26) 



Then 



with 



AM, 



where 
and 



h (n-^(m-)+n]^(m|))-i (m^oE]^(mf)+m^oE]^(m^)) 



A 



A-d 



7e + In 47r 



(27) 



A=0 

(28) 
(29) 



(30) 



7 



5.2. DIAGRAMS CONTRIBUTING 

In a generic way the diagrams contributing are 

p-q p-q p-q 



If \ q q : \ q q ,' \ q 



^, ; 



i J 



These diagrams can be calculated in a straightforward way. For instance 
the W diagram in the ^ = 1 gauge gives 

k=l 



k=l 

where Bq and Bi are the Passarino-Veltman functions, and O'^'^", O'^^'^" are 
coupling matrices. Explicit expressions can be found in [6]. 

5.3. GAUGE INVARIANCE 

When calculating the self-energies the question of gauge invariancc arises. 
We have performed a careful calculation in an arbitrary gauge and 
showed [6] that the result was independent of the gauge parameter ^. 

5.4. THE ONE-LOOP MASS MATRIX 
The one-loop corrected mass matrix is 

M^'^=Mfag + ^M^'' (31) 

where 

Mfag=J^MnJ^ (32) 

Now we diagonalize the 1-loop mass matrix 

Mlgg=M'M^'^M''' (33) 
Then the mass eigenstates are related to the weak basis states by 

Xo —■'^ia Xo 
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with 

J\f^L ^ (35) 

The usual convention in neutrino physics 

J^a = Uak Vk (36) 

is recovered in our notation as 

Uak= ■A/4+fe,4+a (37) 

5.5. SOLAR AND ATMOSPHERIC NEUTRINO PARAMETERS 

Assuming hierarchy in the masses rriu^ and m,^3 and neglecting C/e3 (that 
has to be small) we write the usual two neutrino mixing angle as 

sin2(20i2) = 4C/fiC/22 (38) 

and 

W{2e^^) = AU%(l-U%) (39) 

6. Results 

6.1. MASSES AND MIXINGS 

The BRPV model produces a hierarchical mass spectrum for almost all 
choices of parameters. The largest mass can be estimated by the tree level 
value as shown in Fig. (1). As the figure shows, correct Am^^^ can be 
easily obtained by an appropriate choice of |A|. The mass scale for the solar 
neutrinos is generated a 1-loop level and therefore depends in a complicated 
way in the model parameters. This is shown in Fig. (2) where we have 
fixed the SUSY parameters. The parameter e^/|A| is the most important 
in determining the solar mass scale, but some other parameters also play a 
role [6]. 

Now we turn to the discussion of the mixing angles. We have found that if 
e^/|A| <C 100 then the 1-loop corrections are not larger then the tree level 
results and the flavour composition of the 3rd mass eigenstate is approxi- 
mately given by 

C/«3 ~ A«/|A| (40) 

As the atmospheric and reactor neutrino data tell us that Vt oscilla- 

tions are preferred over z/„ — > i^g, we conclude that 



Ae < A„ ~ A^ 



(41) 
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Figure 1. The atmospheric Am^ as function of \K\/ {^^ M-z^) 




10-2 10"^ 10° 10^ 10^ 10^ 



eV|A| 

Figure 2. Neutrino masses as a function of e^/|A| 

are required for BRPV to fit the the data. This is sown in Fig. (3). 

For the solar angle the situation is more complicated and there are two 
cases to consider [6]. With the usual SUGRA assumptions, ratios of Cj/ej 
fix the ratios of Aj/Aj. Since atmospheric and reactor data tell us that 
Ag -C A^,A,- in this case only the small angle solution can be obtained in 
the BRPV model as shown in Fig. (4 a). However we have shown that even 
a tiny deviation of from universality of the soft parameters at the GUT 
scale relaxes this constraint. In this case the ratio ej/ej is not constrained 
and also large angle solutions can be obtained as shown in Fig. (4 b). 
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6.2. CONSEQUENCES FOR THE ACCELERATORS 

One of the attractives of the BRPV model is that besides accommodating 
the solar and atmospheric neutrino data it can make definite predictions for 
accelerator physics. Then the model can be tested. As R-parity is violated, 
the neutralino is unstable. For this to have experimental consequences the 
neutralino has to decay well inside the detector. This is indeed the case 
as shown in Fig. (5 a). We have seen before that the ratios |Aj/Aj| and 
€i/epsilonj \ were very important in the choice of solutions for the neutrino 
mixing angles. What is exciting is that this ratios can be measured in 
accelerator experiments. In Fig. (5 b) we show the ratio of branching ratios 
for semileptonic LSP decays into muons and taus: BR{x I^q'q) /BR{x 
Tq'q) as function of A^/A,-. We can see that this is directly correlated with 
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Figure 6. Ratio of branching ratios for semileptonic (a) and leptonic (b) LSP decays 
into and t as functions of the BRPV parameters. 

atmospheric angle. In Fig. (6 a) is shown the ratio of branching ratios for 
semileptonic LSP decays into muons and taus: BR{x eq'q)/BR{x 
Hq'q) as function of Ag/A^. This is directly correlated with C/g3. Finally in 
Fig. (6 b) we show the ratio of branching ratios for leptonic LSP decays 
into muons and taus: BR{x — > ver)/ BR{x ^fJ-T) as function of |ee/e^|. 
This is directly correlated with solar angle. 

7. Conclusions 

We have shown that there is a viable model for SBRP that leads to a very 
rich phenomenology, both at laboratory experiments, and at present (LEP) 
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and future (LHC, NLC) accelerators. Most of these phenomenology can be 
described by an effective model with explicit R-Parity violation. We have 
calculated the one-loop corrected masses and mixings for the neutrinos 
in a completely consistent way, including the RG equations and correctly 
minimizing the potential. We have shown that it is possible to get easily 
maximal mixing for the atmospheric neutrinos and both small and large 
angle MSW. We emphasise that the lightest neutralino decays inside the 
detectors, thus leading to a very different phenomenology than the MSSM. 
If the model is to explain solar and atmospheric neutrino problems then 
many signals will arise at future colliders. Thus the model can easily be 
proved wrong. 
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